Statistics of finite-time Lyapunov exponents in the Ulam map 
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The statistical properties of finite-time Lyapunov exponents at the Ulam point of the logistic 
map are investigated. The exact analytical expression for the autocorrelation function of one-step 
Lyapunov exponents is obtained, allowing the calculation of the variance of exponents computed 
over time intervals of length n. The variance anomalously decays as 1/rP . The probability density of 
finite-time exponents noticeably deviates from the Gaussian shape, decaying with exponential tails 
and presenting 2"~^ spikes that narrow and accumulate close to the mean value with increasing 
n. The asymptotic expression for this probability distribution function is derived. It provides an 
adequate smooth approximation to describe numerical histograms built for not too small n, where 
the finiteness of bin size trimmes the sharp peaks. 
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INTRODUCTION 

In a dynamical system, Lyapunov exponents (LEs) 
quantify the average exponential rate of expansion or 
contraction of infinitesimal volume elements in each di- 
rection of phase-space. Negative exponents imply con- 
vergence of initially nearby trajectories, while positive 
ones mean exponential divergence of neighboring orbits 
therefore signaling chaotic behavior. 

By the theorem of Oseledec [1,2], assuming ergodic- 
ity, each exponential rate of growth converges to a time 
asymptotic limit, Aoo, independently of the particular 
trajectory chosen (for almost all initial conditions). A 
useful generalization of asymptotic LEs are the local 
finite-time LEs, A„, calculated over a time interval of 
length n along a given trajectory [3]. In contrast to 
asymptotic ones, exponents defined for finite time de- 
pend on the initial conditions. Local finite-time expo- 
nents must be distinguished from other local exponents 
also considered in the literature such as the so-called local 
finite-sample exponents [4]. The analysis of finite-time 
LEs is particularly important as they are the quantities 
actually measured in computational studies, unavoidably 
performed for finite time. Local LEs have proved to be 
useful in the characterization of a variety of phenomena 
ranging from two-dimensional turbulence [5] to the so- 
called Loschmidt echo [6] . These local exponents are fluc- 
tuating quantities that may even change sign depending 
on the degree of heterogeneity of the relevant subset of 
phase space. Their fluctuations give rise to a nontrivial 
probability density P(A„) that asymptotically collapses 
to a Dirac delta function centered at Aoo- The devia- 
tions from this limit and the convergence to it contain 
rich information on the underlying dynamics [3,7]. 

A relevant question on the dynamics of chaotic sys- 
tems is the existence or not of true trajectories lying 
close to the numerically generated ones. It is known that 
finite-time LEs fluctuating around zero cause the nonex- 
istence of such shadowing trajectories [8]. Particularly, 



in the synchronization transition of lattices of coupled 
chaotic one dimensional (ID) maps [9], unshadowability 
of chaotic orbits is due to unstable dimension variabil- 
ity, characterized by the second largest LE fluctuating 
about zero [10]. The connection between the probabil- 
ity distribution function (PDF) of finite-time LEs in the 
synchronized states and the PDF of the local exponents 
in the uncoupled map [11] motivates the study of the 
statistical properties of finite-time LEs in chaotic maps. 
Within this context, the numerical nature of the difficul- 
ties involved makes analytical results specially relevant. 

This work focuses on the statistical properties of finite- 
time LEs for the logistic map at a parameter value yield- 
ing fully developed chaos, namely, at the Ulam point, 
where x i— > 4x(l — x). This particular mapping allows ex- 
act calculations that are expected to be useful for more 
general systems belonging to the same class of chaotic 
behavior. We begin by calculating the autocorrelation 
function of one-step LEs. This function allows to evalu- 
ate exactly the variance of local LEs. Furthermore, we in- 
vestigate thoroughly the PDF of finite-time LEs. Taking 
as a starting point the work by Prasad and Ramaswamy 
[12], we extend it by deriving an user-friendly expression 
for the PDF of finite-time LEs, a very good approxima- 
tion for not too small n. 



FINITE-TIME LYAPUNOV EXPONENTS 

For a ID map x i— > f{x) with differentiable /, finite- 
time LEs A(n), computed over the interval of length n, 
are given by [7] 
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where Xj — /•'{xq) is the j*'' iterate of the initial con- 
dition xq. We will compute local LEs after a transient 
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has elapsed, so that averages over different reahzations 
can be calculated with the weight given by the invariant 
density p{x) of the chaotic attractor. 

Random initial conditions yield fluctuations in the val- 
ues of A„ around (A„) = Aqo, with variance cr^(n) = 
(A^) — A^, that by means of (1) can be expressed as 



- IL — J./t — J. 
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where 

Cj,, = (ln|/'(x,)| ln|/'(xfe)|)-(ln|/'(x,)|)(ln|/'(x,.)|> 

(3) 

is the two-time autocorrelation function of one-step (uni- 
tary time) LEs [13]. 

Separating diagonal from off-diagonal terms, one gets 
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If stationarity can be attained, then temporal transla- 
tional invariance holds. In such case, correlations depend 
on the two times only through their difference, so that 

n— 1 n— 1 n—1 

EE^^-.fc = E["-"^]c'M' (5) 

j=0 k>j m=l 

where C(m) = Co,m, and additionally C(0) = Co,o = 
Cm,m, Vm. Hence, 

1 9 ""^ 

a^{n) = -C(0) + —Y^[n-m]C{m). (6) 
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Higher moments would require the calculation of corre- 
lation functions involving several times. 

In order to evaluate the variance explicitly, one must 
compute the autocorrelation function C characteristic of 
the particular mapping / considered. In what follows, 
we will perform the explicit calculations for the logistic 
map f{x) = iix{l — x) at /i = 4, with x G [0, 1] (usually 
called Ulam map). 



TIME CORRELATIONS IN THE ULAM MAP 

For this logistic map the dynamics is exactly solvable 
(see, for instance, [3,15]). In fact, using the change of 
variables z = 2 arcsin(\/x)/7r, it is easy to find that 
the j*'^ iterate is xj = sin^(2-'~'^7rzo). The time evo- 
lution for the new variable is given by the tent map 
^ 1 - 2]^; - 1/2|, with z G [0,1] [3,15]. The invari- 
ant density, p{x) = [7r-\/a;(l — .x)]"'", becomes imiform 
for z G [0,1]. In the z-representation, the required cor- 
relations can be straightforwardly obtained by following 



the procedure of Nagashima and Haken [16] to calculate 
the time autocorrelation function Cjf (m), 

Cnim) = {H{g^iz))H*{z))-{H{g"^{z))){H*{z)}, 

(7) 

of an arbitrary function H{z) where z evolves through 
the mapping g(z) and * denotes complex conjugate. The 
basic idea is to consider the Fourier expansion 

oo 

H{g"^{z))= J2 a^e"^''^ m = 0,1,2,..., (8) 

k— — oo 

where 

a^= [ dzH{g'"{z))e-^''''"^. (9) 
Jo 



In our case H{z) = In |4cos(7rz)|, then it is indifferent if 
we evolve z cither with the tent map or with the Bernoulli 
shift g{z) — 2z(mod 1) that will generate the same se- 
quence H{g^"{z)). We will consider the latter case in 
order to use the results of [16]. Recursively, they arrive 
at 
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otherwise. 



Then, using (7), (8), (10) and additionally considering 
that {H{z)) = {H{g"'{z))) = ag, Vm, the correlation 
results 



Now, for our particular case we have 
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Substituting these coefficientes into (11), we obtain the 
following exact expression for the stationary autocorre- 
lation of one-step LEs: 



C(m) 



12 ' 



tt" 1 
24 2'" ' 



if m = 
otherwise. 



(14) 



The autocorrelation function decays exponentially fast 
with time, with a characteristic time that is the in- 
verse of Aoo = In 2. Although the sequence {xj} gen- 
erated by the Ulam map is (5-correlated (more precisely. 
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Cx{m) = (xjXj+rn) - {xj){xj+rn) = 'g'^Om [17]), Corre- 
lations of higher powers of Xj do not decay to zero im- 
mediately [18], leading to the exponential decay of the 
autocorrelation of one-step LEs. 

Once the autocorrelation (14) is known, from Eq. (6), 
we obtain 



a'{n) = 1- — 

V / e„2 V on 
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(15) 



This result valid for all n is in agreement with the 
outcomes of numerical experiments [12-14] as well as 
with an asymptotic expression previously found [14]. 
It is noteworthy that the variance decays with time 
as Then the quantity A„ = n(A„ — Aqo) = 

In ]/'(a;j)/2] does not evolve diffusively, since its 
variance tends asymptotically to a constant value. If the 
terms of the sum in Eq. (1) were independent, then the 
decay of cr^(n) should be as 1/n. However, those terms 
are highly dependent: The sum can be expressed as a 
function of the single variable Xq through the map suc- 
cessive iterates (xq). This is a nice example illustrating 
that linear correlation and dependence are quite different 
statistical concepts. Naturally, although the correlations 
decay exponentially fast, the central limit theorem will 
not hold in this case. 



PDF OF FINITE-TIME EXPONENTS IN THE 
ULAM MAP 

Given a PDF px{x), the standard formula for the PDF 
Priy), such that y = (j>{x), is 



PY{y)= / px{x)5{y - (t){x)). 



(16) 



Prasad and Ramaswamy have already presented [12] 
the exact although hermetic expression for the PDF 
P(A„) of local LEs on the Ulam map, that comes im- 
mediately from the knowledge of the invariant measure 
p{x), using (16). In fact, from Eq.(l), one can write 

A„ = ln[G„(a;)]/n, with Gn{x) = 4"n;>o |1 - VKx)V 
Then, one has 



P{Xn) 



dxp{x)5{Xn-HGn{x)]/n), (17) 
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P{\n) = nexp(nA„) 
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(18) 



where the siimmation runs over all real roots of the 
(2" — l)-order polynomial Gn{x) — exp(nA„). Because 
Gn{x) is odd order, there is always at least one real 
root. For A„ below the mean value, all roots arc real. 
As \n increases above the mean value, the roots leave 
the real axis by pairs. At each point where this occurs, 



G[j(.t) = 0, then a new divergence appears in the PDF 
[12]. The shape of the exact PDF for the scaled expo- 
nent A„ = n(A„ — A(x>) is illustrated in Fig. 1 for different 
values of n. 

Prasad and Ramaswamy proposed a smooth ansatz for 
the discontinuous expression (18), that in terms of the 
scaled variable reads 



P{Xn) 



exp(-[A„[) 
l-exp(-2lA„l) 



(19) 



for A„ G (— oo,nln2]. They based their conjecture on 
the analysis for small n. Here we will show that it is pos- 
sible to advance further, developing expression (18) and 
evaluating it explicitly, therefore, finding analytically a 
smooth expression valid for large enough n. 



1 - (a) n = 2 



1 - (b) n = 5 




1 (c) n = 1 
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FIG. 1. Exact P(A„) as a function of An = n(\n — \oo) 
(full lines), obtained numerically from Eq. (21), for different 
values of n indicated in the figure. As n increases divergences 
become narrower and cannot be appreciated unless one refines 
the mesh increasing the number of data exponentially with n. 
The dashed lines correspond to the smooth approximation 
given by (30). 
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Considering that In |/'(a:;j)| = In | [/"(xq)]'] and 

making the transformation x = sm^{TTz/2), from Eq. (1), 
one gets the compact expression 

A„ = -ln[sin(7r2)] +ln|sin(2"7rz)| = (20) 

where the scaled exponent is expressed as a function of 
the state variable z (from now on the subindex will be 
omitted) evaluated at the instant at which the trajectory 
segment of time length n starts. The first logarithmic 
term is a smooth function in [0, 1] that docs not depend 
on n, while the second term is a periodic function with 
divergences at z = j/2", with j & 2 and < j < 2". 
Increasing n will increase the frequency and the comb 
of divergences will become more compact. The shape of 
Fn{z) is illustrated in Fig. 2, for n = 5. 

Since the invariant density of the attractor as a func- 
tion of z is uniform, using (16), the PDF for A„ is given 

by 



P(A„) 



dz <5(A„ - Fniz)) 



E 

j 



1 



(21) 



where z* are such that Fn{z*) = A„. 



to the sign of A„. 

(1) For Xn < 0, the exact PDF is smooth. There is 
symmetry around z = \. In the interval [0, ^] there are 
2" — 1 intersections at z*^ = i/2" ± 5j such that < 
6j « j 1'^"' ■ Taking into account that the main contri- 
bution to the derivative comes from the periodic com- 
ponent, then \F^Xz*±)\ ^ 7r2"|cos(2"(5j7r)/sin(2"(5j7r)|, 
where, from (20), | sin(2"(5j7r)| ~ exp(A„)| sin(j7r/2")|. 
Substitution of the derivatives in (21) leads to 



4e^ 



7r2" 



E 



sin(j7r/2'^ 



^ ^1-62^" sin2(j7r/2»)| 



(22) 



Replacing the summation by an integral by identifying 
X = j/2", we get 



^ ^0 -^1 - e2^n sin2(7ra;) 



(23) 



that gives 



P(A„)~ ^ln[coth(-A„/2)]. 



(24) 




FIG. 2. Scaled Lyapunov exponent A„, = Fn{z) as a fune- 
tion of z, the value of the state variable when the compu- 
tation of the exponent is initiated (full lines). In this case 
n — 5. The dotted line corresponds to the smooth compo- 
nent — ln[sin(7r2)] that does not depend on n. Increasing n 
by one will double the number of periods. There is symmetry 
around z = 1/2, then only half of the full abscissae interval 
has been represented. 



The approximations here performed are expected to be 
good even for small n, namely n > 2. In fact, a good 
agreement between the approximate and exact PDFs in 
this region is observed already for n = 2 (see Fig. 1). 

(II) The domain < A„ < n In 2 is more tricky. At 
each local maximum of Fn{z), a divergence in the PDF 
appears, as it follows from (21). Then, in this region, the 
PDF possesses 2"~^ spikes. Let us consider the interspike 
intervals (A~, A"*"] determined by two successive values of 
z at which F^ = 0. In each interspike interval the PDF 
monotonously increases, because \F'^\ decreases with in- 
creasing A„ (see Fig. 2), and diverges at the upper limit 
of the interval A+, where F^^ = 0. These behaviors can 
be observed in Fig. 1, for n — 2 and 5. As n increases, al- 
though the cusps increase exponentially in number, they 
get thinner, in such a way that they will become invisi- 
ble in numerical histograms due to the finiteness of bin 
size. Then, if n is large enough, histograms would look 
smooth although the exact PDF is not. unless one could 
make the bin size arbitrarily small. In such case, the 
higher the resolution, the larger the height of the spikes. 
However, in practice, there are constraints to reach an 
arbitrarily fine resolution due to finiteness of numerical 
data sets. This fact turns it useful to find a smooth ex- 
pression to describe the histograms actually observed for 
not too small n. In order to do so, we have to calcu- 
late the probability Jj^.^ dXP{X) associated to a finite bin 
of the order of (or larger than) the interspike distance 
AA = A+-A-. 

For each diffeomorphic branch of F„ one has Fl^{z) = 



Let us separate the analysis into two cases according Pni^n ^i^n)) — l/[ 



dXn 



], then, from Eq. (21), the 
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probability that A„ belongs to the interval (A ,A+] is 
given by 



/ dXP{X)=2j2\F-\~X+)-Fr\x-) =2J2A 



(25) 



where Azi arc the segments displayed in Fig. 3, corre- 
sponding to the inverse images of AA. The summation 
runs up to a value m that depends on the particular in- 
terval chosen and that is given by the smooth component 
ofF„. 




Azj Az. 



FIG. 3. Schematic representation of the scaled Lyapunov 
exponent A„ as a function of z. The segments on the horizon- 
tal axis represent the contributions to the probability associ- 
ated to the interspike interval AA = A"*" — A~ . 



This procedure will give rise to a discrete PDF where 
the integrated probabilities (25) can be attributed, for 
instance, to the midpoints A* of each interval AA, such 
that 



P(A*).4A4^=2f:.^^' 



AA 



AA 



(26) 



Now, ~ \F^{z*)\, where the derivative is evaluated 
at the midpoint z* of the interval Azi. As the derivative 
F^ is dominated by the periodic component in Eq. (20), 
one can perform an analysis analogous to that for case 
(I), obtaining 



PCX*) ^ V 



A* J" 



|sin(j7r/2")| 



l-e2^* sin2(j7r/2")| 



(27) 



where jmax — 2" arcsin[exp(— A*)]/7r, for A* G (A~.A+). 
Replacing the summation by an integral as in case (I), 
one obtains 



4eA* re{x*) 
P{X*) ^— dx 



sin(7ra;) 



y'l - e2^* sin2(7ra;) 



(28) 



where 9{X*) = arcsin[exp(— A* )]/?:. After integration, 
one arrives at 



P(A*) ~ ln[coth(A*/2)] . 



(29) 



Let us recall that, in contrast to case (I), here we have 
a discrete PDF defined for the middle points of the in- 
terspike intervals. However, since the interspike distance 
A A goes to zero for increasing n and not too large A„, 
one can smooth the staircase function (29) by extending 
it to every point of the interval (0,nln2]. In that case, 
one gets the symmetric counterpart of expression (24). 

Expressions (24) and (29) can be gathered into the sin- 
gle expression 



P{Xn) 



■ln(coth|A„/2|) . 



(30) 



This result explains the apparent symmetry of numerical 
PDFs for large n despite the fact that left and righthand 
wings of the distribution have very different genealogies. 
For large |A„|, the PDF decays exponentially, yielding 
straight lines in the semi-log representation displayed in 
Fig. 4. Notice that the approximate distribution (30) for 
the scaled variable A„ does not depend on n. In par- 
ticular, it is expected to be exact in the limit n — > oo. 
That is, (30) is a good representation of the delta func- 
tion to which the finite-time distribution asymptotically 
collapses. Remarkably enough, the approximations that 
led to (30) did not spoil the normalization condition, i.e., 
dXp{x) = 1 holds. 

In Fig. 4 the approximation (30) is compared to an his- 
togram obtained from numerical iterations of the map- 
ping. The analytical prediction is in very good accord 
with numerical results and it describes the numerical 
data more accurately than the ansatz (19). Notably, the 
variance of the scaled exponent A„ — n(A„ — Aqo) calcu- 
lated with our PDF is o-^ = tt^/G. Then for exponent A„ 
one obtains a'^{n) = n"^ /{Qn?) which coincides with the 
exact result (15). if neglecting the term of order 1/2" (the 
order of the performed approximations). In contrast, the 
ansatz of ref. [12] 
0.79 7rV(6n2), i.e.^ 
value. 

As an aside comment, notice that in the exact distri- 
bution the ordinate at the origin increases with n but is 
finite for finite n. However, expression (30) diverges at 
A„ = (as also (19) does). In the exact PDF, divergences 
appear for positive abscissae only and accumulate close 
to the origin as n increases. Therefore, strictly speaking, 
a divergence at the origin is expected in the limit n — > oo 
only. 



gives (j\n) = (7rVl2+ [In2]2)/n2 ~ 
more than 20% smaller than the true 



5 



1.0 
0.5 



0.0 




FIG. 4. P(A„) as a function of A„ for n = 10. Small sym- 
bols correspond to the histogram accumulated over 10'^ histo- 
ries obtained through numerical iterations of the map, after a 
transient of 10^ time steps. The bin size is 0.02. Expressions 
(30) (full line) and (19) (dotted line) are plotted for compari- 
son. Inset: the same data in semi-log representation allow to 
appreciate the tails. The spikes in the righthand tail of the 
histogram are visible for large A„ where they become more 
sparse and are amplified by the logarithmic scale. 



SUMMARY AND FINAL REMARKS 

We have derived analytical results for the statistics 
of finite-size exponents on the logistic map at the Ulam 
point. In particular, wc have found the autocorrelation 
function of one-step LEs, which allows to obtain the ex- 
act analytical expression (15) for the variance of local 
exponents, putting into evidence the origin of its anoma- 
lous time decay. Moreover, we have found the smooth 
analytical PDF (30), expected to be the exact one in the 
limit n — > oo. Therefore, we have extended the work 
initiated by Prasad and Ramaswamy [12], providing a 
complete and consistent picture of the statistical proper- 
ties of finite-time LEs at an outer crisis. 

As the special shape of the distribution of finite-time 
exponents at the Ulam point is also seen in other systems 
presenting fully developed chaos [12], our results are ex- 
pected to be useful for more general systems falling into 
the same class. In addition, these analytical results could 
give insights on features occurring in the synchronization 
transition of lattices of coupled chaotic maps [9], such 



as hipcrbolicity breakdown through unstable dimension 
variability [10], since the distribution of finite-time LEs 
in the synchronized states of such networks is related to 
the distribution of local exponents in the uncoupled map 
[11]. 
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knowledge Brazilian agencies FAFERJ and PRONEX for 
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